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Using the Balitsky-Kovchegov (BK) equation as an explicit example, we show that nonlinear 
QCD evolution leads to an instability in the propagation toward the infrared of the gluon transverse 
momentum distribution, if one starts with a state with an infrared cut-off. This effect takes the 
mathematical form of rapidly moving traveling wave solutions of the BK equation, which we investi- 
gate by numerical simulations. These traveling wave solutions are different from those governing the 
transition to saturation, which propagate towards the ultraviolet. The infrared wave speed, formally 

] infinite for the leading order QCD kernel, is determined by higher order corrections. This mecha- 

nism could play a role in the rapid decrease of the mean free path in the Color Glass Condensate 

■ scenario for heavy ion collisions. 
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I. INTRODUCTION 



The properties of perturbative Quantum Chromodynamics (QCD) at high energy and/or density presently receive 
> much attention. On phenomenological grounds this domain is interesting both for deep inelastic scattering at HERA 
and hard scattering at the Tevatron and LHC, and for heavy-ion collisions with an initially large density of partons, 
. e.g. at RHIC and LHC. On theoretical grounds the high energy and density domain is a challenge, since it implies 
a resummation of logarithms in the perturbative expansion at all orders of the strong coupling constant. This 
resummation goes beyond the one at the origin of the linear evolution equation in rapidity Y known as the Balitsky- 
Fadin-Kuraev-Lipatov (BFKL) equation [1] . It implies non- linear terms which describe the damping of the scattering 
\ amplitudes due to the high density of partons (mainly gluons) produced during the evolution. This leads to a transition 
to saturation [2] where a new state of QCD matter, the "Color Glass Condensate" (CGC) is expected to appear [3]. 
Indeed, the CGC is proposed as the initial state formed in heavy ion collisions at high energy [4]. 

In its mean-field version, the QCD evolution towards saturation is described by the Balitsky-Kovchegov (BK) 
nonlinear equation [5, 6]. In the color dipole picture of QCD [7, 8], it governs the energy dependence of the onium- 
target amplitude, where the target is supposed to be interacting independently with each dipole resulting from the 
evolution of the onium wave function with energy. In parton language, and neglecting impact parameter dependence, 
the BK equation describes the rapidity Y evolution of the "unintegrated gluon distribution" , i.e. the distribution of 
gluon transverse momenta k in the target. 

An interesting set of results on this distribution during the transition to saturation have already been obtained from 
the mathematical properties of the BK equation. Indeed, analytic asymptotic solutions of the non-linear equation 
have been found [9] in terms of traveling waves, i.e. scaling solutions T(k, Y) ~ T (k/Q s (Y)) , where L = log(fc 2 
and Y play the role of space and time, respectively. In fact, the traveling wave pattern of the solutions of the BK 
equation, which is already present at pre-asymptotic energies [10], reflects the earlier proposed [11] "geometric scaling" 
of the 7*p cross section. 

The saturation scale Q S (Y) is directly related to the movement in time of the traveling wave and in fact 
dlagQ s (Y)/dY ~ v (v > 0) is nothing else than the speed of the wave which reaches a universal critical value 
at high energies, depending only on the kernel [9]. Hence, in physical terms the gluon transverse momenta propagate 
with rapidity towards higher and higher values, i.e. to the ultraviolet, the typical transverse momentum scale being 
the saturation scale Q s . This is the key property of the transition to saturation. 

In fact the "final" result of this evolution — the Color Glass Condensate — is supposed to be a fully saturated state 
where all gluons have a momentum equal to (or possibly larger than) the saturation scale characterizing the medium. 
However, in the rapidity evolution from an onium-target scattering, this "final" state is not expected to occur easily 
at finite energy for the whole system. On the other hand, heavy ion collisions at high energies could, due to the 
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FIG. 1: Sketch of traveling wave solutions for the propagation towards the infrared and the ultraviolet. For the ultraviolet 
momentum distribution T, the wave front connecting the large-T and small-T travels from the left to the right as Y increases. 
For the infrared momentum distribution U, the traveling wave, when it exists (see section V), travels in the opposite direction, 
from larger to smaller values of the minimum momentum. The functions T ,U are represented for three increasing rapidity 
values. 

high initial density of partons, be a good candidate for the physical occurrence of the CGC [4] in the early stages 
of the collision. However, it is known that difficulties occur in reconciling this appealing scenario with the rapid 
thcrmalization which seems to occur in the next stages of the process. More precisely the mechanisms by which the 
original CGC, being inherently perturbative, feeds the low momentum region [12] are currently under scrutiny [13]. 
Hence, this asks the question of the propagation towards the infrared, when one starts with an initial state with a lower 
transverse momentum limit, or at least characterized by a typical high transverse momentum, as is the case for the 
CGC. This propagation has been considered from the point of view of a bottom-up mechanism [12], possibly helped by 
plasma instabilities [13]. However, it is worth considering the effect of energy evolution upon this type of dense initial 
phase already in the initial stage of the collision. This is our physical motivation for looking at the properties of the 
BK equation from a new point of view, the BK equation being considered as a prototype for non-linear perturbative 
QCD evolution. We expect most of our results to be valid also for more complete nonlinear evolution equations. In 
the present paper we essentially stick to the relevant mathematical properties of the BK evolution, leaving a more 
detailed analysis of the physics of heavy-ion collisions for the future. 

We will thus discuss how the BK equation gives information on the propagation with energy of gluon momenta 
towards the ultraviolet and the infrared; namely, how the gluon momentum distribution evolves towards larger and 
smaller momenta, respectively, when starting with appropriate initial conditions. We will consider initial conditions 
with a "sharp" cut-off, either ultraviolet or infrared (see Fig.l), but the mathematical universality properties of 
traveling wave solutions will allow much more general initial conditions, provided the initial fronts are steep enough. 

The plan of our paper is the following. In section II we recall the BK equation in transverse momentum space. In 
section III, we discuss the propagation of momenta towards the ultraviolet (i.e. large gluon momenta) and show that 
we recover the same traveling wave solutions and scaling properties as for the solution of the BK equation for the 
physical amplitude. In section IV, using appropriate initial conditions with an infrared cut-off, we derive a new type 
of scaling solutions of the BK equation for the propagation towards the infrared (i.e. towards small gluon momenta). 
In section V we solve the same equations numerically, confirming the analytical predictions and study in detail the 
solutions. Conclusions on the propagation properties of the BK equation and its possible relevance for the CGC and 
the initial stage of the thermalization process are discussed in section VI. 



II. THE BALITSKY— KOVCHEGOV EQUATION 



The BK equation for the onium-target S-matrix in position space reads 



dS_ 

&Y 



(x 01 ,Y)= d 2 x 2 JC(x ,x 1 ;x 2 ) {S(x 02 , Y)S(x 12 , Y) - S(x i, Y)} 



(1) 
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where 



/C(x ,xi;a;2) = - ^ — — ( 2 ) 



"01 X Y1 



is the BFKL kernel [1]. Eq. (1) is the impact-parameter independent form of the BK equation, which is what we 
consider in this work. a s is the QCD coupling constant, Xi are the endpoint coordinates of the QCD dipoles in two- 
dimensional transverse space, and Xij — Xi — Xj. To obtain the equation for the amplitude, perform the replacement 
S -> 1 - TV to get 



d/V 
W 



{xoi,Y) = Jd 2 x 2 JC(x ,x 1 ;x 2 ) {M{x Q2 , Y)+N{x 12 , Y)-Af(x 01 , Y)-M{x 02 , Y)N(x 12 , Y)} . (3) 



Again assuming impact-parameter independence, it is well-known [6] that the Fourier transformation to momentum 
space leads to a rather simple form of the evolution equation which is suitable for the derivation of useful analytical 
properties [9]. Thus one introduces the function 

N(k,Y) = f p^ e - ik - x °^(x 0U Y) = [°° -Jo(fcr) Af(r, Y), (4) 

where r = \xo\\, and we have integrated over the azimuthal angle. Since AT only depends on the magnitude r of 
the vector xoi, N depends only on the gluon transverse momentum magnitude k. N{k,Y) can be related to the 
unintegrated gluon distribution of the target, and thus gives a direct information on the distribution of transverse 
momenta of the gluons probed in the target when the total rapidity is Y. 
The equation (3) then takes the form 

^^ = ^[K^N(k,Y)-N 2 (k,Y)] (5) 

where the kernel K is the integral kernel of the BFKL equation in momentum space. It acts on the gluon distribution 

as 



K®N{k,Y) = f 
Jo 



12 



dk 



k' 2 N(k',Y)- k 2 N(k,Y) k 2 N(k,Y) 



\k 2 -k' 2 \ V4fc' 4 + k 4 



(6) 



A comment useful for the further analysis is in order. It appears easier to discuss the propagation towards the 
ultraviolet and the infrared using formulas (5,6) in terms of gluon momenta. Indeed, this is the formulation for which 
the mathematical results in terms of traveling waves [9] are better known. However, the universality properties of the 
traveling wave regime in momentum space are expected to be valid in position space too. Hence, the propagation of 
momenta towards the ultraviolet (infrared) will correspond to the evolution of QCD dipoles towards smaller (larger) 
sizes. 

The mathematical argument for this statement is, briefly, that it is easy to realize (e.g. by means of a Mcllin 
transform [1, 7]), that the linear part of the equation is symmetric by Fourier transform, interchanging k, k' with 
1/r, l/r' in Eq.(6). It is then known from general mathematical properties [9, 14], that the traveling wave regime 
does not depend of the form of the non-linear term even if modified by Fourier transformation. 



III. PROPAGATION TOWARD THE ULTRAVIOLET 



Let us first consider the problem of the propagation of the distribution of momenta towards the ultraviolet, in an 
"universal way", that is in a way independent of the precise physical structure of the target. In order to formulate 
the problem, we want to consider initial conditions defined only by a sufficiently sharp ultraviolet cut-off. In fact, 
interestingly enough, we will sec that the resulting traveling wave pattern is independent of the initial cut-off profile. 

We start with a simple choice of initial condition, namely a ^-function. As we shall see, the asymptotic solution is 
the same for any other initial distribution with a sharp enough cut-off 1 . Hence we choose 

T(k,Y o )=0(-L) (7) 



We check this point numerically in Section V. 
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where L = log (k 2 /k 2 ) , with fco being the initial ultraviolet cut-off. The notation T is now for the solution of the BK 
equation starting from the cut-off (7). When the rapidity increases, the transverse momenta of the gluon distribution 
evolve towards larger and larger k. As we shall see, the solution takes the form of a traveling wave, the same that 
was found [9] for the physical amplitude N (k, Y) itself. The wave front defines a critical region in k moving forward 
with rapidity. 

Let us sketch the arguments in the present case of an initial ^-function. One introduces the Mellin transform with 
respect to the variable k/k . In the long-range tail of the distribution, where T(Y) ~ 0, the non- linear term in (5) 
will remain small when evolving in Y and can at first be neglected. The solution of the linear part will thus remain 
valid. The solution is written 

T (k,Y) ~ / |I i e-iL+uMY-Yo) ; (8) 



c 



where 



w(7 ) = ^ x ( 7 ) = ^ [2^(1) - ^(7) - - 7)] (9) 



is the Mellin transform of the kernel 2 and the prefactor I/7 comes from the initial 9- function condition (7). We are 
expecting a wave front and thus search for a specific velocity v+ characterizing its forward motion. Hence we define 
a comoving frame following the ultraviolet front and defined by the shift L — > Lwf = L — v + Y. A priori, no special 
velocity v+ is selected, since for each value of 7 one has a velocity v+(j) = which stabilizes the value of Lwf- 
However, it can be shown from general grounds [14] that, once the non-linearities are properly taken into account, 
and provided the initial conditions are steep enough in L, a critical velocity 

, + = ^ = min^^( 7c ) (10) 

7c 7 7 

is selected, in analogy with the group velocity in wave physics [15, 16]. Here 7 C 0.6275 is the critical anomalous 
dimension identical to that found for the BK equation [9, 17]. In fact the only condition [9] on the initial boundary 
dependence is T(k,Y ) < e~ lcL . This can also be inferred from analyticity properties of the integrand of Eq. (8) in 
Mellin space [9] , where the point is to check that the position of the prefactor pole - is situated to the left (in 7) of 
the critical point 7 C = 0.6275. This ensures that the overall behavior will be driven by the critical slope 7 C and not 
by the one corresponding to the initial condition. These conclusions have been tested numerically in Ref. [18]. 

Indeed, the argument is essentially the same as for the transition matrix clement, even if the initial condition is 
different. Recall that in that case, the initial conditions are dictated by the perturbative QCD property of "color 
transparency" T(k,Y) oc k~ 2 at large k (i.e. color transparency, in terms of dipoles). In mathematical terms, all 
solutions verifying a sharp enough cut-off T(k, Y ) < e~ lcL lie in the same "universality class" and thus follow the 
same traveling wave pattern at high enough rapidity. For instance, both the physical amplitude N and the solution T 
with ^-function initial condition lie in this "universality class" . Hence a unifying description of "propagation towards 
the ultraviolet" is correctly defined. We shall discuss the "propagation towards the infrared" within the same context 
in the next section IV. 

There exists a critical value k c (Y)/ko which characterizes the position of the moving wave front T of momenta 
propagating towards the ultraviolet. The value of k c (Y) is defined as the value of k, at a specified Y, where the 
distribution T has a fixed (and, for consistency, small) constant value. Since this is also the region where these 
distributions vary significantly with L, see Fig. 1, k c (Y) represents the most probable domain value for momenta as 
a function of Y. On general grounds, one finds 

k c (Y) cx k cxp (v+Y - ^- log Y^j , (11) 

where 7 C and v+ are given by the solution of the equation (10). The correction <~ \ogY is a universal retardation 
effect due to the non-linearities. There exists also a third 0(Y^ 1 ^ 2 ) sub-asymptotic universal term in the exponential 
[9] which we neglect for our purposes. We also leave for further study the interesting question of pre-asymptotics and 
early scaling [10, 18]. 



2 In the expression (8) of the Mellin-transform one may see that it is symmetric in the Fourier transform from momentum to position 
space, allowing to extend to the latter the traveling wave properties of the former, as quoted in section II. 



■5 



Together with the front velocity dk c (Y)/dY, it is possible to obtain some general knowledge of the front profile in 
the region next to the tail, including information on the diffusive approach to the scaling curve. This means that the 
gluon momenta do not reach instantaneously the scaling behavior of the traveling wave. There is a kind of diffusive 
retardation effect, violating geometric scaling. One finds 



my) ~ const - x v^^ los (^)) (Si) %xp (-^o-^ iog2 (^)) • (i2) 

A physically important remark follows from the identification of the traveling wave characteristics of general ultravi- 
olet evolution with the amplitude Af with perturbative initial conditions. It means that the dynamics are dominated 
by gluon momenta of order of the saturation scale rather independently of the ultraviolet cut-off profile (if sharp 
enough). In other terms, the propagation towards the ultraviolet is essentially determined by the evolution of the 
saturation scale. These gluons are expected to eventually build the CGC phase. 



IV. PROPAGATION TOWARD THE INFRARED 



Let us imagine that at Y — Yq, instead of the boundary condition (7), one starts with boundary conditions for 
another solution U of the BK equation, that describes an initial distribution of gluons with a low momentum cut-off 
fco; see 3 Fig. 1. As in the previous case, we shall discuss an ideal 4 ^-function condition, but our results will be valid 
for the whole class of initial conditions lying in the same "universality class" of sharp enough infrared cut-offs, to be 
defined with precision later on. 

We start with the initial condition given by 

U(k,Y ) = 0(+L), (13) 

where U (k, Y) denotes the corresponding solution of the BK evolution equation. We will now demonstrate that under 
very general conditions such an initial condition leads to a gluon distribution that propagates rapidly "backwards" , 
i.e. towards the infrared. 

However, a special treatment is needed for the BFKL kernel given by (2) in momentum space or (9) in Mellin space. 
The critical value of the speed of the traveling wave is obtained from the equation v — w'(7 c ) = ^{lc)/lc, an equation 
that has one finite and one infinite solution when \ is the leading logarithmic (LL) BFKL kernel. Indeed, \ nas one 
unique finite value of j c is j c = 0.6275, and one formally infinite one corresponding to the pole at 7 = 0, see Fig. 
2(a). The infinite value (formally for 7 = 0) is the indication of an infrared instability of the propagation which is 
to be regularized, as will become clear later on. The analysis of the asymptotic solutions of U requires a particular 
treatment in this case. 



A. Diffusive approximation 

Let us proceed in steps and first go back to the original diffusive approximation of the BK equation where known 
mathematical properties are useful to address our problem. This approximation of the LL BFKL kernel, considered 
in Ref. [9], consists in expanding the kernel around 7 —\ and keeping terms up to second order, 

X(7) - x(§) + X'(l)(7 - \) + ^X"(l)(7 - \f (diffusive approx.). (14) 

This brings [9] the BK equation into the form of the Fisher-Kolmogorov-Petrovsky-Piscounov equation [19]. In this 
kernel there are no poles and the equation for 7 C has two solutions, one (7+) with positive, the other (7-) with 
negative value t>_ of the speed, see Fig. 2(b). The positive value yields the same forward traveling waves propagating 
with speed v + discussed in the previous section. 

We now emphasize the significance of the second, negative, value. It will correspond to "backward" traveling waves 
moving with speed V- towards the infrared. 



3 We keep for convenience the same value fco in, e.g., Fig. 1 as for the initial ultraviolet cut-off, but this has no impact on the results. 

4 In the physical systems we have in mind, there should be also a damping factor at fc — > 00. This is required to give a meaningful 
distribution in position space. Due to the universality properties, these limitations are expected to play no important role for our 
conclusions. 
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FIG. 2: Graphical solutions of the equation x(lc) = JcX'ilc), for (a) the LL BFKL kernel, Eq. (9), (b) the diffusive approxi- 
mation, Eq. (14), and (c) the two-pole approximation, Eq. (16) for 70 = 0.2. For the BFKL kernel there is no solution giving 
a negative velocity (corresponding to negative j c ), while for the diffusive and two-pole approximations there are both positive 
and negative velocity solutions. Note that the magnitude of the negative velocity (i.e. the slope of the tangent) is much larger 
than for the positive velocity (in the LL BFKL case, it is infinite). 

The existence of two types of traveling wave solutions for the same equation has already been noticed in Rcfs. [20]. 
However, in that case they correspond to a statistical physics problem for which both waves propagate in the same 
direction 5 but with different speed. In QCD evolution, the speed v_ = x{l-)/l-i W1 th 7- < and x{l-) > is 
negative. As in the forward case, there is a universality class of solutions with the same pattern, provided the initial 
condition for backward propagation is steep enough in the backward direction, namely U(k, Y ) < e 7 - L ,L < 0. This 
universality class is characteristic of the propagation towards the infrared when using the diffusive approximation of 
the BFKL kernel. 



B. Regularization by energy-momentum conservation 



We now consider cases where the initial LL BFKL kernel gets modified, see Fig. 2(c). This is expected e.g. from 
energy momentum conservation constraints and other next-to-leading logarithmic (NLL) BFKL effects [21, 22]. The 
pole at 7 = gets shifted to a negative value (compare Fig. 2-a and -c). Indeed, the effective kernel function x{l) 
used in phenomenological calculations [22] gets modified along these lines and thus allows for a finite negative value 

7— 

A simple argument for understanding the shifting of the pole is energy momentum conservation: In the relation 
w = (a s A r c /7r)x(7), energy conservation translates into the condition that u) = 1 when 7 = 0. In the limit 7^0, the 
BFKL kernel behaves as x ~ 1/7: so we ma Y incorporate energy conservation by writing 

X(7) -ItTT ( sma11 7), (15) 

7+^ 

which moves the 7 = pole to 7 = —a s N c /n. 



C. Regularization using the two-pole model 

It is well-known that the LL BFKL kernel is well-approximated by a two-pole model, x(l) = 1/7+1/(1 — 7); which 
reproduces the poles of the BFKL kernel. This kernel has also been studied in Rcf. [18] in the context of saturation 
and is found to give results very similar to those from the full kernel. 

In view of the above we use as a toy kernel 

x{l) = —^— + -, ^^+4 flog 2--^— ), (16) 

7 + 7o 1-7 + 70 V 1 + 2 7o / 



5 We interpret this difference between the statistical physics problem and QCD by the fact that the former is a one-dimensional problem 
analogous to the random breaking of a segment into smaller pieces [20], while the latter is related to the cascading of dipolcs in a 
two-dimensional position space (even if the equation under study is its one-dimensional projection). 
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FIG. 3: Propagation of momenta for the LL BFKL kernel. Left: towards the infrared; Right: towards the ultraviolet. The 
rapidity steps are SY — 1. In the left plot, all lines with Y / scale on the same curve. 



which reproduces the main features of the BFKL kernel with the poles at 7 = 1 and 7 = shifted by an amount 70 and 
—70 respectively. We consider 70 as a free parameter of our model. This kernel can be viewed as an approximation 
of some next-to-leading effects, and is also simple enough to use in our reasoning. It will allow for analytical and 
numerical 6 studies of the infrared propagation for given 70 and its instability when 70 — > 0. The constant term is 
chosen to reproduce the value of the full BFKL kernel at 7 = 1/2. 

For example, choosing the value 70 = 0.2 as illustrated in Fig. 2(c), one obtains two real solutions to the equation 
defining j c , 7+ = 0.766 and 7" = —0.103, leading to traveling wave solutions with both positive and negative 
velocities. All these examples, including the 70 = case where the infrared instability appears, will now be discussed 
by comparison with numerical simulations. 



V. NUMERICAL SIMULATIONS 



The BK equation has been solved numerically by many groups [23] by direct solution of the intcgro-diffcrcntial 
equation. Recently, in Ref. [18], it was solved using a discretization of the integrand through Chebyshev approximation 
and a Runge-Kutta method for the ensuing system of ordinary differential equation (see [18] for more details on the 
method). The solutions were studied in detail and compared to analytical results obtained using the traveling wave 
front method. In this paper we use the same numerical method to study the solutions for the distributions T and U 
in fc-space. All results that follow are obtained with a fixed strong coupling constant, a s N c /Tr = 0.2. 

First, let us consider the simulation of equation (5) made using the original BFKL kernel (6), corresponding to the 
case of Fig. 2(a). 

We observe that traveling wave solutions are generated in the forward direction, while the backward direction 
is characterized by a fixed (non-traveling) scaling curve, obtained almost immediately in the rapidity evolution. 
Moreover, the solution extends also immediately towards very negative values of L = log(/c 2 /fcQ), i.e. in the very far 
infrared. This is a clear sign of infrared instability of the evolution with the LL kernel. 

In LL BFKL, hence, only forward traveling waves are present. The backward solution propagates instantly to high 
Y, and generates a fixed scaling solution. In a formal way, the speed gets infinite and the wave front is situated at 
L = -co. This situation is illustrated in Fig. 3. 

In order to analyze in more detail this infrared instability, we considered the same problem with much smaller steps 
in rapidity, see Fig. 4. We see the evidence for a rapid diffusive evolution towards the scaling curve. Note also the 
far infrared region, which is immediately reached by the evolution. This is a feature related to the singularity of the 
canonical BFKL kernel at 7 = 0, and, as we shall see now, is physically tamed by next-to-leading corrections. This 
sensitivity to next-to-leading QCD corrections is a specific and interesting feature of backward propagation towards 
the infrared. 



See Ref. [18] for the numerical implementation of the two-pole model. 
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FIG. 4: First steps of BFKL evolution towards the infrared. The rapidity steps are SY = 0.1. 
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FIG. 5: Propagation of momenta for the two-pole model. Left: towards the infrared; Right: towards the ultraviolet. The 
curves are drawn for values 70 = 0, 0.2, 0.4 and the rapidity steps are 5Y = 1. 

In Fig. 5, we show the similar plot for the non-linear evolution using the two-pole kernel (16) with different values 
of the parameter 70 which, in our toy model, schematically characterize the amount of next-leading corrections. 
Interestingly enough, the toy model results for 70 = reproduce the same features as for the BFKL kernel, with 
a scaling curve in the backward propagation. For non-zero 70 = 0.2, 0.4, left-moving traveling waves appear, in 
agreement with the analytical predictions. In all cases the speed u_ is larger than for the forward propagation (i.e. 
larger than the ordinary forward BK traveling waves). Hence the infrared instability is made finite but the propagation 
towards the infrared is still much more rapid. In terms of QCD dipoles, one may say that dipoles with large sizes are 
easily produced during the energy evolution, in the "universality class" of backward propagation. The comparison 
between the forward and backward evolutions is exemplified in Fig. 6. 

To round up the numerical results we give an example of the (in) dependence of the traveling wave regime from 
the precise shape of the cut-offs in the initial conditions provided they are sharp enough, either in the forward or 
the backward direction respectively, to satisfy the mathematical requirements for creating traveling waves. We thus 
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FIG. 6: Comparison of forward and backward propagation for 70 = 0.4. The rapidity steps are SY = 1 
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FIG. 7: Propagation of momenta for different initial conditions. Left: two-pole model towards the infrared; Right: BFKL 
kernel towards the ultraviolet. Continuous lines: initial step functions. Dotted lines: initial smoother cut-off (see text). The 
curves are drawn for the BFKL kernel and the rapidity steps are SY = 1. 

compare initial step functions with an initial form [4] of the cut-off profile which is smoother 7 than a ^-function, while 
satisfying the sharpness condition. As typical examples (see Fig. 7) we considered the backward waves obtained for 
the toy model and the forward waves for the original BFKL kernel. 

Finally, we confirm that the speed of the front propagation is correctly given by the analytical formula for asymp- 
totically large rapidities, i.e., we compare the rapidity dependence of the saturation scale of the forward and backward 
traveling waves obtained from the numerical solution to the analytical formula for asymptotically large rapidities. 

The saturation scale k c (Y) is defined as the momentum, at a certain rapidity, where the amplitude reaches a 
particular value k, i.e., T(k c (Y), Y) = k. The overall normalization therefore depends on the chosen value of k and it 
is most convenient to study the logarithmic derivative of k c (Y), where the normalization cancels. Because the shape 
of the front is not completely fixed as it propagates in rapidity, the functional form of d log k c (Y)/dY will also depend 
on k, but this dependence is a subleading effect in Y and vanishes as Y becomes large 8 . Here we choose n = 0.01. 



In practice we used the numerical Fourier transform from r-space of [4] 

l-cxp(-r 2 Q 2 (Y)/41og[e + l/(r 2 A 2 )]) . (17) 

See [18] for examples and a thorough discussion. 
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FIG. 8: Comparison of the analytical form for the derivative in (18) with the results from the numerical simulation, for backward 
evolution (left) and forward evolution (right). The equation analyzed here is given by the the two-pole toy kernel (16) with 
7o = 0.4. Note the order of magnitude difference between the backward (left) and forward (right) propagation speeds. 



The analytical result given by (11) is 

9 log k c (Y) 



V +-T^T7, (18) 



dY ^ 2 lc Y 

where subleasing OiY ~ 3 / 2 ) terms have been omitted for sake of simplicity. Fig. 8 shows the results of the comparison 
between the simulations and expression (18) for the two-pole toy kernel (16) with 70 = 0.4 (the constants are given 
by Eq. (10): v + = 0.7372, u_ = -6.112, 7^ = 0.9157, and 7! = -0.2202). The results indicate that the two results 
tend to the same asymptotic values v± of the derivative, while there are still appreciable subasymptotic effects [18]. 
A study of subasymptotic effects and early scaling, which we postpone for future work, could take benefit from the 
recent work in Ref. [10]. 



VI. CONCLUSIONS AND PHYSICAL IMPLICATIONS 



Let us summarize the results obtained in our study of the gluon momentum propagation (or equivalently, dipolc 
size evolution) during nonlinear QCD evolution at high density /energy: 

(i) The propagation of gluon momenta towards the ultraviolet and infrared with energy correspond to two different 
"universality classes" of solutions of the non-linear Balitsky-Kovchegov evolution equation. They are selected by the 
existence of an ultraviolet (resp. infrared) cut-off in the initial conditions but largely independent of the cut-off profile, 
if sharp enough. 

(ii) The propagation towards the ultraviolet is characterized by asymptotic traveling wave solutions which arc 
independent of the (sharp enough) initial distribution. In particular, they are the same as those exhibited by the 
onium-target amplitude with initial conditions defined by the QCD transparency property. Hence, the rate of evolution 
towards the ultraviolet is universally governed by the saturation scale, generalizing the behavior of the BK dipole- 
target amplitude. 

(iii) The propagation of sizes towards the infrared shows new and distinctive features which have not been discussed, 
to our knowledge, in the literature previously. In the case of the LL BFKL kernel, i.e. in the leading log approximation 
of perturbative QCD, one obtains a scaling distribution (after a short diffusive interlude) which extends almost im- 
mediately towards the far infrared region, revealing an infrared instability. This phenomenon is due to the singularity 
of the BFKL kernel at 7 = 0. 

(iv) As soon as this singularity is shifted or smoothed out, new traveling wave solutions appear, with a (large) speed 
of backward propagation towards the infrared. In particular, next-to-leading corrections to the BFKL kernel should 
give rise to such a phenomenon and determine the value of the critical speed, as tested in different examples. 

(v) The analytical predictions of the traveling wave pattern in both directions (ultraviolet and infrared) have been 
verified by numerical solutions of the evolution equations. 
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Our findings of the new traveling wave pattern towards the infrared within a purely perturbativc evolution framework 
could have sizable physical implications. An example we have in mind is the problem of thermalization of a Color 
Glass Condensate formed in the first stages of a heavy-ion collision. 

Indeed, it is well known that it is not easy to reconcile the formation of a perturbative CGC (at weak coupling) as 
an initial compound state in a heavy ion collision, with the hydrodynamic features of the final state which probably 
require a fast thermalization or, at least, a rapid decrease of the mean free path in the medium. We think that 
a discussion of this mechanism would be interesting for the discussion of thermalization of an initial Color Glass 
Condensate phase during a heavy ion collision process. 

As a further outlook, it will also be worth examining the effect of the growth of the coupling constant during the 
propagation towards the infrared. Our present results have been obtained for simplicity assuming a fixed value of the 
coupling constant. They will have a natural extension to running coupling, using the method of Ref. [9] for solutions of 
the BK equation with running coupling. We expect qualitatively similar results but the quantitative analysis deserves 
to be performed. 

These essentially geometric features may also be modified by a better description of the interaction with the target. 
It is nowadays emphasized [18, 24] that fluctuations could modify the saturation solutions, violate geometric scaling 
or even exert a back-reaction onto the linear regime. What we obtain from our analysis of the "universality class" 
of propagation towards the ultraviolet will probably survive the implication of fluctuations. Even more interesting is 
the investigation of fluctuation effects on the infrared propagation. This will be the subject of forthcoming studies. 
Properties borrowed from statistical physics could help in attacking this problem with new tools. 
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